Jahn- Teller Distortions and the Supershell Effect in Metal Nanowires 
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A stability analysis of metal nanowires shows that a Jahn- Teller deformation breaking cylindrical 
symmetry can be energetically favorable, leading to stable nanowires with elliptic cross sections. 
The sequence of stable cylindrical and elliptical nanowires allows for a consistent interpretation 
of experimental conductance histograms for alkali metals, including both the electronic shell and 
supershell structures. It is predicted that for gold, elliptical nanowires are even more likely to 
form since their eccentricity is smaller than for alkali metals. The existence of certain metastable 
"superdeformed" nanowires is also predicted. 
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The effects of shell filling on the abundance spectra of 
metal clusters have been known for years and explain the 
existence of clusters with "magic numbers," correspond- 
ing to full electronic shells, which are observed more fre- 
quently than others 0. More recently, electronic shell 
and supershell structures in alkali metal nanowires have 
been reported by the Leiden group 0, 0, The conduc- 
tance G was recorded during the breaking of nanowires in 
a mechanically controllable break junction (MCBJ), and 
histograms were built out of thousands of measurements. 
When the temperature is large enough to allow the wire 
to explore the phase space of possible shapes, the shell 
structure shows up as peaks in the histogram of VG that 
are equally spaced, while the supershell structure mani- 
fests itself as a modulation of the amplitude of the peaks 
(See Agrait et al. 5] for a recent review). 

The electron-shell structure in alkali nanowires can be 
understood qualitatively using the free-electron model of 
a cylinder |2| . However, more detailed stability analyses 
of axisymmetric wires plZlla] have revealed a sequence of 
stable "magic" radii, with characteristic gaps that are not 
fully consistent with the nearly perfect periodicity of the 
experimentally observed peak positions. The deviations 
can be accounted for neither by the inclusion of disor- 
der 9||, nor by the use of more elaborate self-consistent 
jellium models |l0|. Since gaps in the sequence of cylin- 
drical nanowires arise from a degeneracy of conductance 
channels, it is natural to assume that a Jahn- Teller defor- 
mation breaking the symmetry can lead to more stable 
deformed configurations. 

While the experimental manifestations of electron-shell 
structure are similar in metal clusters and nanowires, 
the Jahn- Teller effect plays out quite differently in these 
two systems. Metal clusters and nanowires differ funda- 
mentally in that surface effects tend to stabilize clusters, 
while they lead to a Rayleigh instability in nanowires . 
As a result, Jahn- Teller deformations of clusters are very 
common and typically rather small Il2l | while, as we 
show here, they only occur for a minority of nanowires 
and can be rather large. Nevertheless, there are several 



stable nanowires with elliptic cross-sections (referred to 
as elliptical nanowires henceforth), leading to a sequence 
of stable cylindrical and elliptical nanowires that allows 
for a satisfying interpretation of the experiments on elec- 
tronic shell and supershell structure 0,y,0|. In addition, 
our theory predicts the existence of short-lived superde- 
formed wires. 

We use a nanoscale free-electron model, treating the 
electrons as a non-interacting Fermi gas confined within 
the wire by hard- wall boundary conditions fl3j . This con- 
tinuum model is especially suitable for alkali metals, but 
can also be applied to other monovalent metals. Macro- 
scopic arguments Q suggest that the surface tension a s 
should drive a crossover from a crystalline solid to a plas- 
tic or fluid below a critical radius i? m i n = o- s /o~y, where 
cry is the yield strength. Typically, i? m i n ~ 10 nm in met- 
als pfl. This implies that electronic effects should dom- 
inate over atomistic effects for sufficiently small radii. 
Indeed, a crossover from atomic-shell to electron-shell ef- 
fects with decreasing radius is observed in both metal 
clusters an d nanowires 0, 0, 0, 0, justifying a 
posteriori the use of the nanoscale free-electron model in 
the later regime. 

Guided b y e vidence of the existence of non-spherical 
clusters [ill [l2| , we consider straight elliptical nanowires 
[lij aligned along the z-axis. The wire cross-section is 
characterized by the ellipse's two major semi-axes a and 
b, or equivalently, by the area parameter a 2 = a ■ b and 
the aspect ratio e = a/b. When a small perturbation, 
written as a Fourier series, is added to a wire of length 
L, the geometry is characterized by the two functions 

cr(z) = d + \^<j q e iqz = a + XSo-(z), 
i 

e{z) = e + \^2e q e iqz = e + X5s(z), (1) 
i 

where the dimensionless parameter A sets the size of the 
perturbation. The perturbation wave vector q must be 
an integer multiple of 2tt / L and, since o~(z) and e(z) are 
real, we have o~- q = a* and e_ ? = e*. The deformation 



must fulfill another constraint which comes from the fact 
that, depending on material parameters, the deformed 
wire tries to find a compromise between a volume con- 
serving deformation and one ensuring electroneutrality 
[Tij. Here we use the general constraint 

TV = kpV - 7](3Trk 2 F A/8) = const, (2) 

where V is the volume of the wire, A its surface area, and 
kp = 2tt/Xf the Fermi wavevector. The parameter rj can 
be adjusted so as to fix the value of the effective surface 
tension to the material-specific value. In particular, r\ = 

corresponds to a constant- volume constraint, and rj 

1 is the constraint of constant (zero temperature) Weyl 
charge. 

A nanowire is an open system, with electrons being 
injected and absorbed by the leads, so that we need to 
calculate its grand canonical potential Q in order to de- 
termine the energetic cost of a deformation. A fully quan- 
tum mechanical stability analysis of cylindrical nanowires 
H shows that, if the wires are short enough, and/or the 
temperature is not too low, their stability is essentially 
determined by their response to long wavelength pertur- 
bations. For these we can use the adiabatic approxima- 
tion and decouple the transverse and longitudinal mo- 
tions yielding a series of effectively one-dimensional prob- 
lems, which can be solved with the WKB approximation, 
in which f2 is given by |l3| 
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Here / = (1 + exp[(E - E F )/k B T]) 1 is the Fermi func- 
tion at temperature T, Ep is the Fermi energy, and the 
E v 's are the transverse eigenenergies. The sum in Eq. (0} 
runs over all open channels is, for which E v {a, e) < E. A 
factor of 2 accounts for spin degeneracy. Note that at zero 
temperature Eq. Q reduces to Cl[T=0,L] = E[Ep,L]. 

The transverse eigenenergies can be written as 
E v {a,e)=E F 

Jkp a )i i where the 7^'s are the zeros of the 
modified Mathieu functions [2oj| . and can be computed 
numerically. In the limit e — ► 1, we recover the result for 
a cylindrical wire, where most eigenenergies are degen- 
erate. As shown in Fig. ^ this degeneracy is lifted for 
e ^ 1 as a result of the decrease in symmetry, so that an 
elliptic nanowire can sometimes be more stable than the 
axisymmetric one, leading to a Jahn- Teller deformation. 

The energetic cost of a small deformation of a straight 
elliptical nanowire can be calculated by expanding 
Eq. |J2J) as a series in the parameter A, 

Sl = Sl<® + Afi« + x 2 n& +0(\ 3 ). (5) 

A nanowire (a, e) is energetically stable at temperature 
T if fiWfoe) = and n^(a,e) > for every possible 
deformation (Se,5a) satisfying dAf = 0. 




FIG. 1: Eigenenergies of a 2d electron gas confined to an 
elliptical shape as a function of the aspect ratio e. Solid and 
dotted lines indicate even and odd states respectively, while 
states non-degenerate at e = 1 are shown with a dashed line. 



The expansion of Eq. yields 




where E v = E v (a,s), and 7„ = 7i/(e)- The prime indi- 
cates differentiation with respect to e. The elements of 
the matrix A in Eq. (JSJ are given by 
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The constraint dAf — on allowed deformations restricts 
the number of independent Fourier coefficients in Eq. JQl . 
Hence we can express cro in terms of the other coefficients 
and expand it as a series in A, 



^ 0) (e, a,e ) + Xa^ (e,a,{e q ,a q }) +0(X 2 ), (10) 



and eliminate <To from Eqs. (J7J| and |JSJ, modifying the 
expressions for and the stability matrix A. At zero 
temperature, the modified condition = determines 
which wires are stationary states, while the positivity of 
the modified matrix A reveals the stability of the wire. 
The results for finite temperature can be derived analo- 
gously by evaluating Eq. @ numerically. 
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FIG. 2: (color) Energetically stable 
cylindrical and elliptical nanowires as a 
function of temperature (in units of the 
Fermi temperature Tf). The surface ten- 
sion was adjusted to a value of 0.22 N/m, 
corresponding to Na pl|. The aspect ra- 
tio e is coded via the scale shown in the 
inset. The numbered arrows label the 
most stable configurations that enter our 
analysis of the shell structure. Arrows A 
and B label two highly deformed ellipti- 
cal wires that are addressed specifically 
in the text. 



Figure [21 shows the stable geometries as a function 
of temperature and conductance, with the surface ten- 
sion adjusted to the value 0.22 N/m (i.e. ?y = 0.93), 
corresponding to Na |2lJ. The color represents the 
value of the aspect ratio e and the x-axis is given by 
the square root of the corrected Sharvin conductance 
G s = G {k 2 F a 2 /4 - k F V/4ir + 1/6). Here G = 2e 2 /h 
is the conductance quantum and V the ellipse perimeter. 

The possibility to observe a stable wire in conductance 
histograms depends mostly on two conditions: (i) The 
wire has to be formed often enough to be statistically 
relevant; (ii) It needs to have a long enough lifetime so 
as to be recorded. Condition (i) depends on the ability 
of the system to probe new configurations, which is de- 
termined by the mobility of the atoms and the density 
of stable geometries in configuration space. Regarding 
point (ii), a calculation of the lifetime of a nanowire is 
beyond the scope of this paper, but one can expect that 
wires that are linearly stable up to larger temperatures 
T m ax will have longer lifetimes at low temperatures |22| . 

Based on these considerations, we extract the most sta- 
ble configurations from Fig. [2] defined as the geometries 
that persist up to the highest temperature compared to 
their neighboring configurations. Those wires are marked 
with numbered arrows in Fig. [21 For each stability peak, 
we extract its mean Sharvin conductance and its width, 
and plot them as a function of the peak number in Fig. [3J 
together with the experimental data from Ref. 2J. Note 
that the striking fit is only possible when including ellip- 
tical nanowires, for which the corresponding aspect ratios 
e are also shown in Fig. [3 

The heights of the dominant stability peaks in Fig. [21 
exhibit a periodic modulation, with minima occuring at 
peaks 4, 8, 12, 18, etc. The positions of these minima are 
in perfect agreement with the observed supershell struc- 
ture in conductance histograms of alkali metal nanowires 
Interestingly, the nodes of the supershell structure, 
where the shell effect for a cylinder is suppressed, are 
precisely where the most stable elliptical nanowires are 
predicted to occur. Thus the Jahn- Teller distortions and 
the supershell effect are inextricably linked. 

While both the shell and supershell Q effects are 
accurately described by our stability analysis, our ther- 



modynamic model does not directly address the complex 
dynamical process by which various contact structures 
form in MCBJ experiments, and therefore cannot de- 
scribe every detail of the conductance histograms, such 
as their highly non-trivial temperature dependence. 

Our stability analysis also reveals two highly deformed 
elliptical nanowires with conductance values of 2 Go and 
5 Go, marked in Fig. [21 by arrows A and B, respectively. 
They are expected to appear more rarely due to their 
reduced stability relative to the neighboring peaks, and 
their large aspect ratio e (see Table that renders them 
rather isolated in configuration space 23]. Conductance 
histograms of the alkali metals often do show a conduc- 
tance peak at 5 Go and a shallow shoulder at 2 Go Q>9- 
Nevertheless, these peaks were not included in the exper- 
imental analysis of shell structure in Refs. Q and 0. Note 
that the definition of an experimental "peak" depends on 
the selection criteria and smoothing of the data. 

Of the three alkali metals, Li, Na and K, Potassium has 
the lowest and Lithium the highest melting temperature. 
This suggests that at a given temperature the mobility 
of the atoms is highest in K and lowest in Li. This is 
reflected in the fact that the evidence of the highly de- 
formed nanowires A and B is clear for K, can be seen for 
Na, but is barely visible for Li 0. 
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FIG. 3: (color online) Comparison of the experimental shell 
structure for Na (taken from Ref. 0) with our theoretical 
predictions of the most stable Na nanowires. Elliptical wires 
are labeled with the corresponding aspect ratio e. 
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TABLE I: Most stable elliptical nanowires (with e > 1) for Na 
and Au. Listed are the labels from Fig.|5J the mean Sharvin 
conductance Gs, the aspect ratio e, and the maximum tem- 
perature T max up to which the wires remain stable, which is 
related to the depth of the corresponding energetic minimum. 



In addition to wire A, Fig. [5] shows a number of "su- 
perdeformed" nanowires with e > 1.5. The most stable 
are at (G, e) = (2, 1.65), (3, 2.25), (4, 2.8), (6, 1.65), (8, 
2.0), (12, 1.65), (19, 1.55), (29, 1.55) and (41, 1.6). An 
aspect ratio near 2 is favorable for the shell effect |24| . 
as evidenced by the large gaps in the energy spectrum 
(Fig. [I}, but the large cost in surface energy renders these 
wires less stable than the magic cylindrical wires and the 
"normal-deformed" wires listed in Table [I] 

Conductance histograms have also been recorded for 
gold, both at low |2^| and room temperatures (l5lllH . lr^ |. 
where evidence of electronic shell effects has been re- 
ported. Although the explanation of some features of Au 
wires, e.g. surface reconstruction, requires explicit inclu- 
sion of 5d-orbitals, the jellium model is sufficient to ex- 
plain electron-shell effects. Performing a stability analy- 
sis for gold requires adjusting the deformation constraint 
|2l according to the surface tension for Au, (1.3 N/m 
[2l|, T) = 0.61), which is higher than for Na. Therefore 
the results are somewhat different: Stable configurations 
appear at the same conductance values, but the tem- 
perature T max up to which elliptical wires remain sta- 
ble is larger, and these wires are less deformed than for 
Na (see Table QJ. As a consequence, the probability to 
observe elliptical nanowires in experiments is enhanced. 
Indeed, several experimental histograms [lil Esj ] for Au 
show clear peaks at conductance values of 2 Go and 5 Go- 

In conclusion, we have presented a stability analysis 
of elliptical metal nanowires, using a jellium model, and 
have shown that Jahn- Teller-distorted wires can be sta- 
ble. The derived sequence of stable cylindrical and ellipti- 
cal geometries explains the experimentally observed shell 
and supershell structures for alkali metals. Deformed 
wires can explain additional conductance peaks observed 
in alkali metals and gold. 
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